Topological Non— connectivity Threshold in long-range spin systems 
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We demonstrate the existence of a topological disconnection threshold, recently found in Ref. Q], 
for generic 1 — d anisotropic Heisenberg models interacting with an inter-particle potential R~ a 
when < a < 1 (here R is the distance among spins). We also show that if a is greater than the 
embedding dimension d then the ratio between the disconnected energy region and the total energy 
region goes to zero when the number of spins becomes very large. On the other hand, numerical 
simulations in d — 2, 3 for the long-range case a < d support the conclusion that such a ratio remains 
finite for large N values. The disconnection threshold can thus be thought as a distinctive property 
of anisotropic long-range interacting systems. 

PACS numbers: 05.50.+q, 75.10.Hk, 75.10.Pq 



I. INTRODUCTION 

Despite the wide use in statistical physics, long-range 
interacting systems, that is those systems characterized 
by a pairwise interaction decaying as a power law of the 
mutual distance with an exponent a less than the em- 
bedding dimension, do not have a well defined thermo- 
dynamic limit . Also is it not at all clear whether their 
equilibrium properties can be described by the ordinary 
tools of statistical mechanics. For instance, the nonequiv- 
alence between the microcanonical and the canonical ap- 
proach has been recently found in a long-range rotators 
model in the thermodynamic limit 0- 

Besides these relevant implications in the foundation of 
statistical mechanics and in theoretical physics as wellQ, 
the non-extensive behavior of long range systems has 
nowadays become important for applications too, rang- 
ing from neural systems to spin glasses 0. 

Within the class of long-range interacting systems, 
classical spin models, widely investigated during the last 
years 0, are the most easy-to-handle both from the an- 
alytical and the numerical point of view. Within such 
class of systems, (to be more precise, a class of anisotropic 
Heisenberg models) the existence of a threshold of discon- 
nection in the energy surface has been demonstrated^ 
for an interparticle interaction with infinite range. It 
has been called non-ergodicity threshold for historical 
reasons^, even if the term can generate some confusion. 
Indeed non-ergodicity is only an obvious consequence: 
it simply means that the energy surface is topologically 
disconnected in two regions characterized by positive and 
negative magnetization. In other words it cannot exist 
a dynamical path connecting them and all trajectories 
starting from one region of the phase space stay there 
forever. For this reason we prefer here to call it Topolog- 
ical Non-connectivity Threshold (TNT). 

The presence of the TNT cannot be considered an ex- 
otic mathematical peculiarity of some toy model. Its dy- 
namical relevance has been studied in [Jj], where an ex- 
plicit expression for the reversal times of the magnetiza- 



tion (the time necessary to jump from one branch to the 
other) has been given in the neighbors of the critical en- 
ergy point. Reversal times diverge at the TNT as a power 
law with an exponent dependent on the number of the 
particles (and, probably, on the embedding dimension) as 
in ordinary phase transitions. Strictly speaking, even if 
in different context and for different models, the relation- 
ship between energy thresholds and topology transitions 
in the configuration space of classical spin models has 
been recently investigated^^- 

Also, while the threshold was explicitly found within 
a class of anisotropic classical Heisenberg models with 
an easy axis of magnetization and all-to-all constant in- 
teraction, at the same time systems with nearest neigh- 
bor interaction were found to have a different behavior. 
For instance, the portion of disconnected energy region 
grows with the number of particles N, less than the en- 
ergy itself, thus resulting in a zero ratio in the thermo- 
dynamic limit. Needless to say, such ratio stays finite for 
anisotropic coupling and all to all interaction. 

While this feature is surely due to the anisotropy of 
the coupling (such finite ratio disappears for isotropic 
coupling even in the case of infinite range interaction), 
the question arises whether the presence of the TNT 
can be considered a pathological effect of the unphys- 
ical infinite interaction range or it is just somehow re- 
lated with the long-range effects. This does not repre- 
sent an academic question. Indeed, despite the possible 
applications of such model even in the case of all-to-all 
interaction |l l| . physical models require taking into ac- 
count more realistic interactions, usually anisotropic^^ 
and depending generically from the inter-spin distance, 
as for the dipole-dipole coupling or when the spin is cou- 
pled with the electron spin of the conduction band of 
a metal, e.g. the RKKY model 13] . This leads quite 
naturally to Hamiltonians with an interparticle potential 
decaying as a generic power law with an exponent a of 
the relative distance R. The results found in Ref.Q can 
thus be recovered by letting respectively a — > (all to all 
coupling) or a — > oo (nearest neighbor coupling). 
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Here, we extend the previous results to the whole 
class of models with an inverse power distance poten- 
tial and show that a = 1 is a critical exponent for non- 
connectivity in d — 1 chains. 

In general, the extension to higher dimensions is far 
from trivial, both numerically and analytically. However, 
we prove that the TNT, if any, can not "survive" (and 
we will specify the precise meaning below) when N — > 
oo, and a > d. Numerical simulations in 2-d and 3- 
d also suggest that, for a < d, the ratio between the 
disconnected energy region and the total energy range is 
finite in the thermodynamic limit. We thus conjecture 
that the TNT is a generic property of anisotropic long- 
range systems in any dimension. 



A dynamical consequence of the TNT is that below it, 
a sample with a given initial magnetization m y , cannot 
change the sign of m y for any time, since the constant 
energy surface is disconnected in a positive and a negative 
magnetization regions, thus no continuous dynamics can 
bring an isolated system from one region to the other. 

Our proof will follow two steps: in the first part we find 
the minima of the x and y parts separately. Then we will 
show that the disconnected ratio goes to zero for long- 
range interaction, while it goes to some finite constant in 
the short range case. 

III. ONE DIMENSIONAL CASE 



II. THE MODEL 

The Hamiltonian is a simple generalization of that con- 
sidered in Ref. yj, and it is given by: 



1 N 



(1) 



where Si — (Sf , Sf , Sf) is the spin vector with contin- 
uous components and modulus 1, TV is the number of 
spins, rj (0 < r) < 1) is an anisotropic coefficient, and 
c |i-j| = K — il~ Q ' w hh a > 0. For definiteness we con- 
sider here only the case of an even number N of classical 
spins. 

Such kind of models are characterized by a minimal 
and maximal energy E min , E max , and by a finite energy 
range E max — E m i n that we call energy spectrum (ES). 
In order to define properly the disconnection threshold, 
let us introduce the set A of all spin configurations with 
a zero projection of the total magnetization along the 
y-axis: 



JV 



A = {C(5i, . . . Sjv)| m y = 



The TNT is thus defined as: 



0}. 



Edis = Mince a[H], 



(2) 



(3) 



and the spin configurations corresponding to E c u s will be 
indicated as Cdis- Here we are mainly interested in all 
those cases where the TNT, if any, occupies a significant 
portion of the ES in the thermodynamic limit. For this 
reason let us define the disconnection ratio: 



Edis — En 

\E m 

1:11 



> 



(4) 



A system will be considered disconnected only if r — ► 
const. > 0, when TV — > oo. Note that the definition of r 
given in Eq. has a meaning only for systems with a 
bounded energy range. 



A. TNT, if any, is in the XY plane 

Roughly speaking, since Hamiltonian is indepen- 
dent of the z-component of spins, the minimum will oc- 
cur when the spins are as large as possible in Q , namely 
in the XY plane. 

In order to prove that the configuration Cdis effectively 
lies in the XY plane, let us assume that it has some S z 
component different from zero. For definiteness assume 
Sf > 0. It is then possible to define another configuration 
C simply making a rotation around the y-axis clockwise 
or counterclockwise which puts the spin Sf onto the plane 
XY. The energy difference between these two configura- 
tions can be computed at glance: 



AE 



>1 



N 

E< 

i=2 



c i - 1 S?[±Jl-(S%)*-S 



(5) 



Here the different sign ± indicates the different way 
(clockwise or counterclockwise) of rotation. Since Sf = 
iyl — (Sf) 2 — (Sf) 2 , it is then clear that, according to 
this sign it is always possible to rotate in such a way to 
have AE < 0. 

The same procedure can be applied n times for any 
other Sf ^ 0, so that we will end with a configuration 
G A (the rotation does not change the constraint) 
with energy E^ < Edis ■ We can therefore consider the 
configurations in the XY plane. This choice has the main 
advantage that it is sufficient to consider as independent 
variables the angles 9i of the i-th spin w.r.t. the £-axis, 
thus satisfying automatically the conditions on the unit 
spin modulus : Sf — cos 9i and Sf — sin#,:. Therefore 
we have to minimize the following expression: 

H = — ^""^ c \i-j\ (v cos Qj cosOj—sinOisinOj) = rjH x +H y , 

(6) 

under the constraint X)j=i sm @i = 0- Since 

E dis = Min(H\m y = 0) > Min(nH x )+Min(H y \m y = 0), 

_(7) 

a lower bound of Ed is can be provided finding the minima 
of the two terms in the r.h.s. of Eq. (7J. Note that 



the first term on the right of Eq.Q does not contain 
constraints, indeed we will show in the next section that 
the absolute minimum of H x automatically satisfied the 
constraint m v = 0. 



B. Minimum of H x 

The minimum of H x due to the overall plus sign (77 > 
0), can be obtained as in a standard anti- ferromagnetic 
spin system with neighbors interaction, that is disposing 
alternatively the spins along the x-axis as +1 and — 1. 

Indeed, for any a, let us call the k — th spin component 
Sf = s and rewrite the energy as follows: 

E * = s ^2 c \j-k\ s j + 2 c \k-i\ S k S i = as + & > 

where a, b are constants independent of s. There are 
two possibilities, a — or a 7^ 0. In the first case the 
energy E x turns out to be independent of the k — th 
spin, while for o / 0, the minimum is attained when s 
has its maximal value (+1 if a < 0, —1 if a > 0). So, 
in any case we can say that the minimum occurs when 
\s\ = 1. Since the procedure can be iterated for all spins 
components, the minimum occurs when S% = ±1 for any 
k, namely in the class of Ising models (ci, cr/v) with 
<r, — ±1 and long range interaction. 

We still have to prove that the minimal energy is ob- 
tained when the spins have alternating signs. To this end, 
let us consider the interaction between two neighbor spin 
pairs, the j-th, namely <J 2 j-\o 2 j and the (J + k + l)-th, 
namely <72j+2k+xO'2j+2k+2 (this can be done since there 
is an even number of spins). There are 16 possibilities 
but only six of them have different energy due to the 
symmetry on the whole change of sign. They are: 
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= 2 + c 2 fc+i + 2c 2 fc +2 + 


C2fc+3 


+ + 


+- E 2 


= C2/C+1 — C2/C+3 




+ + 


-+ £3 


= -E 2 




+ + 


— £4 


= 2 — c 2 k+i — 2c 2k+2 — 


C2fc+3 


+ - 


+ - £5 


= — 2 — C 2 fc+1 + 2c2fc+2 


— C 2 fc+3 


+ - 


-+ E 6 


= -2 + c 2fe+ i - 2c 2 fc+2 


+ C 2 fc +3 



From the monotonicity of the function c x one gets, 
E5 < Ei,E 2 ,E 3 , E4. 
Imposing E5 < Eq we obtain 

2C2fc+2 < C2fc+1 + C2fc+3, 

that holds from the convexity property of c x . 

Since £5 is minimal for each pair interaction, the ab- 
solute minimum will be obtained using, for each pair, the 
configuration 
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FIG. 1: (Color online) Spin values along the chain vs the spin 
index (only the central part of the chain has been shown) for 
the Hamiltonian H y . a) for fixed a = 0.1 and different JV 
values as indicated in the legend; b) for fixed N = 20 and 
different a values (see the legend). 

The energy E x for this configuration can be computed 
immediately: 

^=EW^. (9) 
fe=i 

Let us notice that such a result is far from being ob- 
vious. Indeed a decreasing non-convex function c x could 
give rise to a different minimal configuration. 

C. Minimum of H y 

Let us now switch to the more difficult task (due to 
the constraint) of computing E y = Min(H v \m y = 0). 
Physically, due to the overall minus sign in front of H y , 
one can expect that clusters of aligned spin with unit 
modulus (ferromagnetic order) will decrease the energy 
with respect to other configurations. This is surely true 
for nearest neighbor interaction (a = 00) but it can not 
be true for all a values. For instance, when a = 0, the 
energy corresponding to the configuration with half spins 
equal to 1 and half equal to —1 (the order is irrelevant) 
is Eq = N/2 > 0, while the true minimum E = is 
attained when all spins are 0. 

Then, the question arises of what can be the minimum 
in presence of a generic a. 

Applying the standard Lagrange multipliers formalism, 
one has to minimize the function: 

H = — c \i-j\ sin6>; sin#j — A sin 0j. (10) 

where A is the Lagrange multiplier associated to m y = 0. 

Taking the derivatives, we get, for each spin, two pos- 
sible solutions: 



c x = {s? = (-i) i }U 



(8) 



cos 0i = 0, 



(11) 



E 



sinfl, — A = 0. 



(12) 



However, solving, even numerically, the system H11I12JI is 
more difficult than finding directly the minimum. 

We have therefore calculated the minimal configura- 
tion under constraint, using an iterative optimization ap- 
proach based on the FFSQP solver and also devel- 
oping the following approach outlined here below. 

• Start with a random configuration with m y = 0. 

• Chose for the k — th spin a new value between —1 to 
1 and compute the energy. This generally produces 
a change in magnetization Am y ^ 0. 

• Distribute equally Am y among the other spins 
taking into account the constraint about their mod- 
ulus. Specifically subtract /add to every spin the 
minimum of its distance from the values ±1 and 
the mean of Am y . 

• Iterate over all spins up to an energy variation less 



than some fixed value (from 10 
simulations). 



to 10 in our 



The two approaches give the same result: for any ini- 
tial random configuration the algorithm described above 
converges to some smooth configuration, for any finite 
N and a > 0, as indicated in Fig. ^ There, we consid- 
ered respectively the case of a fixed varying N (Fig.[IJi), 
and N fixed varying a (Fig. ^>) . Within the numerical 
errors the spins in the minimal energy configuration are 
distributed monotonically and anti-symmetrically w.r.t. 
the center of the chain. Then we assume that E y is given 
by an anti-symmetric distribution of the spin with a non- 
decreasing (or non increasing) monotonic dependence of 
the y— spin component along the chain. 

An interesting feature is the presence of a finite domain 
wall (defined by those spins having length less than 1) 
between two clusters with = +1 (|) and Sf = —1 
(|) respectively. With decreasing range of interaction 
(increasing a: Fig. [Q)) or increasing number of spins 
(Fig. QJi), the interface region between the clusters (f) 
and (|) decreases. This agrees, at least qualitatively, 
with the results obtained for the nearest neighbor model 
(a = oo) where the minimal configuration is given by 

c n = a ... ra - i)- 

Thus, due to long range interaction, an interface region 
between the two clusters with opposite magnetization is 
produced. It is, of course, physically relevant to under- 
stand if the size of the interface region goes to zero in the 
N — > oo limit. 

Strictly speaking the configuration E^ is not an abso- 
lute minimum for any a > and finite N. In order to 
prove that, consider the configuration 




FIG. 2: Critical N cr as a function of a. The region above the 
line is where the TNT is given by E s (one spin pair decreased) 
and the region below it is where the TNT has more than one 
spin pair decreased (E s t). 



with N — 2 spins satisfying condition (|llfl and the two 
central ones satisfying the condition l|12|) . The energy E s 
corresponding to C s can be written as 

„2 



E s = E + cis - 2s[ci - c N / 2 ], 



(13) 



where E is independent of s. The minimum is thus ob- 
tained when s = 1 — c N / 2 / c i 1- The energy difference 
to E n is 



AE = E, 



-N/2 



< 0. 



(14) 



Therefore, for any finite chain and finite a, E s < E^. 

Physically, C s has an energy less than due to border 
effects. Indeed the energy of two opposite spins of length 
1 is E = ci, while it is only c\s 2 for two shorter spins 
\s\ < 1. On the other hand the interaction between s (— s) 
and the other spins of length 1 (—1) is canceled one to 
one but the interaction with the closest spin (— s c\) and 
with the last opposite one (s cat/ 2 ). This gives Eq. lfT^ . 

The same procedure can be applied taking a trial con- 
figuration with energy E st : 



(T ... V,t;s;-s;-t;i ... |). 



In this case a minimum with s < t < 1 can be found 
only for N < N cr (a) = (2 a+1 - 1/2) 1 / 01 . Asymptotically, 
for large N, this implies that for N > 2C 1 / Q , where 
C = 2e~ 1 / 4 > 1, the minimal solution has energy E s . In 
Fig.|21we show the graph of N cr (a), and the two regions 
in the plane (N,a), where E s is the minimal solution, 
and where another minimal solution with four (or more) 
spins with length less than 1 is possible (E st ). Since 
N cr (a) — > oo for a — > 0, for any a/0a sufficiently large 
./V > N cr (a) value exists (thus in the thermodynamic 
limit) such that E s is the minimal solution. Then, for 
N > N cr (a): 



(T...T;s;-s;|...|) ! 



E y - E M 



N 



2a 



(15) 



5 




FIG. 3: (Color online) E x (full black line), E n (dashed red 
line) and Edis (open circles) versus a, for r\ = 0.9 and different 
N values : (a) N = 10, (b) N = fOO. 



where E^ can be written in closed form as: 



. 2\ Q_1 N ^3k-N N/2-k 

+ ~^- + (N/2 + k)°- (16) 



We thus proved that for a > and N > N cr (a) 

E n > E dis > V E X + E y , (17) 

where the expressions for E-^i, E Xl and E y are given 
respectively by Eqs. ifTol 151 fT5|l . 

D. Thermodynamic limit 

Let us now show that, in the long range case < a < 1, 
the ratio r between the disconnected ratio, defined by 
Eq. , goes to a non zero constant when the number of 
spins goes to infinity, while, for short range interaction 
a > 1, it goes to zero, thus revealing the intrinsic long 
range nature of the TNT. 

The minimum energy, having as a configuration 
C m in — {S % y — l}iLx (all spins aligned along the y- 
direction) can be easily found: 



E„ 



N ^k-N 

2-~i ha 



k=l 



(18) 



Let us also define the quantities: 

E y + rjE x - E n 



ri 



r-j, 



|-^mm | 



\E„ 



(19) 



Due to (jnj, < ri < r < r 2 . 



1. Long Range 

Consider first the long-range case < a < 1. The fol- 
lowing asymptotic expression, for N — > oo, can be found 
by substituting sum with integrals: 



E,, 



N 



2-a 



E 



11 



N 



(2 - a)(l - a) 
1 - 2" 



E. r 



(l-a)(2-a) 
baN + OiN 1 -"), 



0(N), 
0(N), 



(20) 

(21) 
(22) 



where b a > is a constant independent of N. 

Since both r x ->• 1 2 - 2 a | and r 2 1 2 - 2 a | for TV oo, 
it follows r — * |2— 2 Q | too, so that the disconnected energy 
region remains finite w.r.t. the ES in the thermodynamic 
limit. This prove the disconnection of the system below 
the TNT. It is interesting to note that, as a — ► 1, r — > 0. 



2. Short Range 

In the short-range case, a > 1, one can write the fol- 
lowing asymptotic expression (by substituting sums with 
integrals) : 



E n 



c a N + 0{N 



2-a\ 



where -1 + 1/(1 - a) < c a < -1 + 2 1 " Q /(1 - a). 
Let us first show that, as N — > oo: 



lim 



Tl 



= 0. 



Computing explicitly the l.h.s. of 124(1 one gets 

N-l 



(23) 



(24) 



2 Z^ 72 " 1 1 
< lim — V r 

y fe=i 



E 

fc=JV/2+l 



-JV/2 



< lim — 

W^oo N 



dx ~ l ~ a 



N 
N/2+1 



, N-X\ 

dx = 0. 



(25) 



Then, r 2 — » and, since r 2 > r > it follows r — * 
and the system is not disconnected. 

This concludes our proof, which is valid for anti- 
symmetric coupling (positive rf). 



E. Id, numerical solution for the full model 

Despite the proof of the existence of the TNT did not 
require the explicit knowledge of the spin configuration, 
it may be of some interest to find it. 

Finding analytically the spin configuration of the full 
model under the constraint m y = for any a, rj and N is 





FIG. 4: (Color online) (a) x and y spin components for the 
numerical TNT. Here is N = 100, 77 = 0.9 and a = 0.05. 
(b) x and y spin components for the numerical TNT. Here is 
N = 10, 77 = 0.9 and a = 00. 



a complicated task. Indeed, depending on the different 
values of the parameters, the minimal configuration can 
completely change its shape, for instance from all spins 
along the x-axis with alternating signs (giving rise to the 
energy E x ) to all spins along the y-axis (first half positive, 
second half negative) giving rise to E^. For instance, 
when a = 0, E x < < JS-fj., while for a — > 00 and N 
sufficiently large E^ < E x < (for small N it is also 
possible to have E x < E^ < 0). 

This is explicitly shown in Fig. [31 where E x , E^ and 
Efiis (the TNT) obtained numerically for two different N 
values have been plotted as a function of a. As one can 
see, for a less than some value depending on 77 and N, say 
ao(N,rj), one has E x < E^, while for a > a (N,r]), one 
can have different possible situations depending on the 
77 and N values. For instance for N — 10 and 77 = 0.9, 
E x < En (Fig. Ek) for a — ► 00, while E n < E x for 
a -> 00 and N = 100 (Fig. EJd). 

It is also instructive to describe the behavior of Edi s 
as a function of a. As one can see (Fig. |3Ji), for rela- 
tively small a, Edi S closely follows E x while for large a 
values, Edis, even if different from both, is closer to E^ 
than to E x . This is the rule, at least for large N values, 
and the difference between the configurations given by 
E^i and Edi S is only restricted to a small domain wall 
in the central part of the chain, see Fig. where the 
configuration Cdis in a long range case has been shown. 
Completely different is the situation for small N values, 
e.g. Fig. 0Jd. Here only 10 particles are considered. In 
this case, for a — > 00 (see discussion above) E x < E^i 
and the configuration Cdis is something between C x and 
C n (see Fig. Eb.) _ 

Let us analyze in detail the behavior of the domain 
wall under a change in the parameters of the system. To 
this end we define the energy domain as 



FIG. 5: (Color online) Domain wall energy as a function of JV 
for different a values, as indicated in the legend, and rj = 0.9. 
Also shown as horizontal dashed line |A_E'| = rj 2 /{l + rj). 



Its behavior for different N and a values has been shown 
in Fig. O As one can see, in the large N limit, the energy 
domain approaches, for any a, some finite nonzero value 
dependent only of 77. This is remarkably different from 
the domain wall obtained by minimizing H y under the 
constraint m y = 0, see Sec. IIII CI where the formation of 
the domain wall was essentially due to border effects and 
whose energy disappears for large N values, see Eq. <|15ll . 

This asymptotic value can be understood as follows: 
consider the trial configuration 



r — 

^xy — 



{o,. 



.0. • \ T 

,1,+s,- 



s, — 1. 



-vT 



,0,. 



,0} 
(27) 



II ( ............. - ........ 1} . 

The energy E xy {s) of this configuration is given by 

E*y(s) =E + Cl s 2 - 2s[ Cl - c N/2 ] - 77Ci(l - s 2 ), (28) 

where E is independent of s and rj. 
The minimum, as a function of s, is 



ci + 



(29) 



Then, for N — > 00, s m i n — > 1/(1 + 77) , which is inde- 
pendent from a. This value can be compared with our 
numerical results. The energy difference to in the 
limit N — » 00 is: 



AE' = E xv (s m i n ) — E 



n 



>r 



1 



v 



< 0. 



(30) 



E,, 



\Ed 



E 



1 ! 



(26) 



Its absolute value has been indicated as a horizontal 
dashed line in Fig. [S] As one can see, all curves are 
close to | AS' | even at N ~ 100. While we cannot exclude 
that other configurations, with four or more central spins 
Sf < 1, have an energy less than E xv (s m i n ), we surely 
have found a minimal configuration whose energy is dif- 
fering from E^i for a finite quantity in the limit N — > 00. 
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2D 



3D 



FIG. 6: (Color online) TNT for the nearest neighbor interac- 
tion (a = oo) vs the parameter r\. Here is N = 20. Different 
energies are indicated in the caption inside the figure. 



This in turn implies that the domain does not disappear 
in the thermodynamic limit. 

As a last remark, let us stress that the no-disconnection 
of a system does not necessarily imply the non existence 
of a TNT different from the minimal energy at finite N, 
even for short range interaction. To this end we con- 
sider the strongest short range coupling, namely the near- 
est neighbor one (a = oo) and compute numerically the 
TNT. Results are presented in Fig. (JJ. 

As one can see the numerically computed Edis is differ- 
ent from E m i n for rj ^ 1 so that a finite range of energies 
E m in < E < Edis exists for finite N and nearest neighbor 
interaction. Increasing N, the size of this energy range 
remains constant, while E m i n ~ N. That is why the ra- 
tio r — > for large N values. From the same Fig. El it 
is clear that Cdis goes continuously from a configuration 
close to Cfj. when r] <C I to one close to C x when r] ~ 1. 



IV. MULTIDIMENSIONAL CASE 

The results obtained in the previous Sections for d = 1 
can be extended in greater dimension d > 2. 

While it can be easily shown that, in the short-range 
case d < a, the system cannot be disconnected in the 
thermodynamic limit, the proof of the disconnection 
for the long-range case is essentially based on the as- 
sumption that the minimum energy with the constraint 
m y — is given by an obvious extension of what we have 
found in d = 1. This is, at the moment, justified only by 
our numerical simulations. 

Let us then consider a d-dimensional hypercube of side 
L, such as L d = N and divide it in two equal halves. Let 
us then put half of the spins with y-component in one 
region and the other half in the remaining with opposite 
y-component and call E^ the resulting energy for such 
configuration. Surely the TNT has an energy value less 
or equal to Efi, that is Edis < ^Ti- 
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FIG. 7: (Color online) Critical qo as a function of the number 
of spins for different 77 indicated by different colors, (a) d = 2, 
(b) d = 3. 



Let us then write: 



Emin — El + El + Vf\ , 

E U = E] + Ei + V<[i, 



(31) 



where E^ , E± are the energies of the respective halves 
and V\i, Vf| are the interaction energies between the two 
halves with, respectively, antiparallcl and parallel spins. 

Since E-\ = Ei and — Vj-f = V\i > 0, one has: 



< r < r n = 2 



(32) 



We will make use of the results found in Ref . [TjJ , that 
in our variables read as: 



lim E ^(d a ,N) = 

for a 7^ d and where the constant Cd(a) > for d < a 
and Cd{a) < for d > a depends only on d and a. 



1. Short-Range 

Let us discuss the short range case a > d. In this case 
we have, 



E min (d,a,N) 
N lun oo N = Q(a) ' 

and, since E-\ — E m i n (d, a, N/2), we can write 



(34) 



n<- <- , 2E^/N-E mm /N 

° " " ~ rn = 2 \E mm /N\ ►O.farJV-.oo. 

(35) 

This proves that, in the short-range case, r — > for 
N -> 00. 



8 



/ 
\ 
t / 
\ \ 
/ j 
N 1 




1000 



FIG. 8: (Color online) TNT for the 2-d square lattice, for the 
case a — 0.1, r\ — 0.5, TV = 16 x 16 = 256 spins. Here is 
E dis = -655.968997, while E n = -654.7308. 



FIG. 9: (Color online) Domain energy as a function of the 
number of lattice spins, for the 2-d square lattice and 77 = 0.5, 
asterisks (a = 0.1), circles (a = 1). For a > 1, Edomain 
becomes smaller than the computer precision. 



2. Long-Range 

In the long range-case, a < d, let us assume that, for 
large N values, E dis — > E^. 

Estimate <|33[) becomes in this case: 



,. E min (d,a,N) 
hm — — = C d (a), 



AT-kxj pp-a/d 

so that, for N — > 00, 



(36) 



rn = 2 



2E^/N 2 - a / d - E mm /N 2 - a / d 



2 - 2 a l d . 



\E mm /N 2 -«/ d \ 

That way, r — > const ^ for N — > 00 and a ^ d, and 
a finite disconnected energy range exists in the thermo- 
dynamic limit. 

It is also interesting to note that, as a — ► d, r — ► 0, so 
that the result (|35|l is recovered. 

The disconnection of the system in the long range case 
can thus be proved if we assume Cdis ~ Cu- 

Numerical simulations confirm this assumption. In- 
deed, let us define ceo(r),N) as the smallest value such 
that, 

E x (a 0) rj,N) = E n (a ,n,N). 

Its general dependence on parameters has been presented 
in Fig.[7|for d = 2, 3. As one can see, cto ~ 1/N — > when 
— > 00. In the same picture we indicate the regions 
where E x or E^ are respectively the minimal energies 
satisfying the constraint m y — 0. Also, cto = is plotted 
as a dashed horizontal line, showing that the short-range 
case (above the line) is characterized by E^, while the 
long-range case (below the line) can have different behav- 
iors (E x or E^i), even if physically interesting long-range 
interactions are generally characterized by E^. 



As for the spin configuration Cdis, both for d = 2 and 
d = 3, all physically significant cases can be represented 
by C-fj_. Deviations occur for a << d, where a domain 
wall appears, see for instance Figs. 0D As one can see, 
Cdis is generically represented by two macroscopic blocks, 
with opposite sign of the y-magnetization, with a do- 
main wall at their interface. In the domain wall the y- 
components increase in absolute value toward the center 
and the x-components are more or less distributed with 
alternating signs, see Figs. |H1 

Defining the domain energy, as the difference between 
the numerically found TNT and the energy (J2HJ), one 
has that with increasing N it goes to some constant or 
zero value (see Fig. so that, in the thermodynamic 
limit Edis ~ -E'Tf) wn i cn justify, at least numerically, our 
previous assumption. This concludes our proof of the 
disconnection of the system with 77 > for long-range 
interaction in any dimension d. 



V. NEGATIVE 77 

In this last part we briefly discuss the case 77 < 0. First 
of all 77 > is not a necessary condition for the existence 
of a finite disconnection region. 

Let us first consider the 1-d case; In Eq. @ r]H x be- 
comes ferromagnetic as H y , and the configuration 



C x = {s x = \}f =1 e A, 



(38) 



has an energy E' x < E x (which is the energy of the con- 
figuration C x , see Eq. JSJ)). Moreover, since the number 
of parallel spins in C' x is larger than in Cj j,, we will expect 
that the energy E' x , even if decreased by a factor 77, will 
become soon or later less than E^. 

While this has no consequences in the case a > 1, (we 
still have r-i — > and the TNT does not exist), in the 
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FIG. 10: -Edis (symbols) as a function of 77 for the long range 
case a = 0.5 and AT = 64 (d = 1), AT = 8 x 8 (d = 2), 
iV = 4 x 4 x 4 (d = 3). Also shown as dotted horizontal lines 
En, and, as dashed transverse lines E' x . 

long-range case (0 < a < 1) some interesting features 
appear. 

From Eq. O one has n -> |2 + ry - 2 a |/2 for — > 00 
and a finite disconnection energy region still exists for 
2 Q — 2 < 77 < 0. 

In the other case — 1 < 77 < 2 — 2" nothing can be said, 
even if, according to our numerical simulations Cdis ~ 
This effectively happens in all dimensions d — 1,2,3, as 
indicated in Fig. I1UI where E<n s as a function of 77 as 
been shown in a long range case. As one can see Edi S is 
close to E' x (equal, within numerical accuracy, according 
to our simulations) for 77 < r) cr (a,N), while it becomes 
close to En, for 77 > r) cr (a,N). That holds true in any 
dimensions. 

Let us note that, as a realistic ?7 cr , we can assume the 
intersection point between E' x and En, see Fig. (|1(JII . An 
estimate, that holds only in the thermodynamic limit can 
be obtained following the considerations made in Sec. (|IV|) 
giving r] cr ~ 1 - 2^ . 

It is also clear that assuming: 

Edis ^]E X f]E m i n 



one has r — > 1 + 77, and the system is disconnected even 
for negative 77 in all dimensions. 



VI. CONCLUSIONS 

Summarizing, we have studied the occurrence of 
a topological non-connectivity threshold (TNT) in 
anisotropic Heisenberg models in d = 1,2,3 with an inter- 
action strength depending on a power law of their relative 
distance with the exponent a. We have found that the 
system, in the thermodynamic limit, is disconnected only 
in presence of a long-range interaction < a < d. On the 
other side, in the short-range case, the ratio between the 
disconnected energy region and the total energy region 
goes to zero when N — > 00. The anisotropy represents in 
this class of systems a necessary condition : indeed, in 
the isotropic case, the TNT coincides with the minimal 
energy, thus there is no disconnected energy region. 

Future investigations concern the experimental evi- 
dence of TNT, for instance by looking for the divergence 
of de-magnetization times[j| as a function of energy in 
small magnetic samples. 

Finally, let us point out that from a quantum mechan- 
ical point of view the classical disconnection does not ex- 
clude the nipping of the magnetization through Macro- 
scopic Quantum Tunneling [16J |. Thus the existence of 
TNT could give the possibility to study the emergence 
of Macroscopic Quantum Phenomena in a wide energy 
range (for macroscopic long range interacting systems), 
as has been shown in 0] , where the quantum signatures 
of the TNT in an anisotropic Heisenberg model with all- 
to-all interaction have been studied, and the relevance of 
the TNT w.r.t. Macroscopic Quantum Phenomena ad- 
dressed. 
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